GAMMA-RAY BURSTS AS THE DEATH THROES OF 
MASSIVE BINARY STARS 
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ABSTRACT 



It is proposed that gamma-ray bursts are created in the mergers of double neutron 
star binaries and black hole neutron star binaries at cosmological distances. Bursts with 
complex profiles and relatively long durations are the result of magnetic flares gener- 
ated by the Parker instability in a post-merger differentially-rotating disk. Some bursts 
may also be produced through neutrino-antineutrino annihilation into electrons and 
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■^J- ■ positrons. In both cases, an optically thick fireball of size < 100 km is initially created 
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which expands ultrarelativistically to large radii before radiating. Several previous ob- 
jections to the cosmological merger model are eliminated. It is predicted that 7-ray 
bursts will be accompanied by a burst of gravitational radiation from the spiraling-in 
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binary which could be detected by LIGO. 



Subject Headings: Accretion — Black Holes — Gamma Rays: Bursts — Gravitation 
- Magnetic Fields — Neutrinos — Pulsars — Stars: Binaries — Stars: Neutron 
— Stars: Supernovae — X-rays: Binaries 

Submitted to Ap J (Lett.), March 24, 1992 

1 Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 

2 Princeton University Observatory, Princeton, NJ 08544 

3 Racah Institute for Physics, The Hebrew University, Jerusalem, Israel 



1. INTRODUCTION 



Recent results obtained with the Burst and Transient Source Experiment (BATSE) 
on the Compton Gamma-Ray Observatory (Meegan et al. 1992) suggest strongly that 
gamma-ray bursts (see Higdon & Lingenfelter 1990 for a review) originate at cosmo- 
logical distances. The 153 bursts reported so far appear to be isotropic in the sky and 
to have a distribution of V/V ma , x that is consistent with a cosmological population 
extending to redshifts z ~ 1 (Mao & Paczyhski 1992, Piran 1992). The required event 
rate is ~ 10~ 6 yr -1 per L* galaxy, and the typical energy released in a burst is ~ 10 51 
erg. The short rise times observed imply a source size ~ 100 km. 

Based on these requirements, many authors (Paczyhski 1986, Goodman 1986, 
Eichler et al. 1989, Piran 1990, Narayan, Piran & Shemi 1991, Paczyhski 1991, Piran, 
Narayan & Shemi 1992) have suggested that 7-ray bursts arise in the merger of binaries 
consisting of either two neutron stars (NS-NS) or a black hole and a neutron star (BH- 
NS). This suggestion is attractive for a number of reasons. (1) The model employs 
a known source population; four NS-NS radio pulars have been discovered in the 
Galaxy, viz. PSR 1534+12, PSR 1913+16, PSR 2127+11C, and PSR 2303+46. (2) 
The scenario invokes orbital decay through the emission of gravitational radiation, for 
which we have direct observational evidence in the case of PSR1913+16 (Taylor & 
Weisberg 1989); moreover, three of the above four pulsars have merger times shorter 
than the Hubble time. (3) In a merger, an energy > 10 53 erg will be released in a time 
~ 1 ms and within a radius < 100 km, satisfying the observational constraints. (4) 
NS-NS and BH-NS mergers are estimated to occur at the rate of ~ 10~ 6 — 10~ 5 yr -1 
per galaxy (Narayan, Piran & Shemi 1991, Phinney 1991), in good agreement with 
the observed burst rate. 
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Many arguments have been made against cosmological scenarios in general and 
the merger model in particular. Among these, one objection appears at first sight to 
be quite serious. If 10 51 erg of 7-rays are created in a volume of size 100 km, the 
optical depth due to 7 + 7 — > e + + e~ will be extremely large and the photons will 
apparently be trapped (Schmidt 1978). This objection was refuted by Paczyhski (1986) 
and Goodman (1986) who showed that an optically thick ball of energy, a "fireball," 
will expand relativistically and thereby radiate most of its energy when it becomes 
optically thin. Because of relativistic beaming, a distant observer receives a burst of 
radiation whose temperature and duration will be similar to the initial temperature 
and initial light-crossing time of the fireball. Relativistic beaming also circumvents 
the so-called "Ruderman limit" (Ruderman 1975), which sets an upper limit to the 
distance of a source for a given source temperature, flux, and variability timescale. 
(Beaming solves a related problem in the case of extragalactic radio jets). 

Although the problems associated with the large optical depth and the Ruderman 
limit have been solved, several other objections remain. (1) How is the energy converted 
into 7-rays? (2) How can one avoid baryon contamination which will significantly 
modify the evolution of the relativistic fireball? (3) How can one obtain bursts with a 
median duration ~ 10 s when the dynamical timescale is < 1 ms? (4) Why are burst 
profiles so complex and individually unique? (5) Why do some bursts have a precursor 
several seconds before the primary burst (Murakami et al. 1991)? (6) What produces 
the power-law 7-ray spectrum, which extends well beyond the pair creation limit of 
511 keV and sometimes even up to a few hundred MeV? (7) How can one explain 
the cyclotron absorption features seen in a few bursts? (8) What is the explanation 
for the redshifted electron-positron annihilation lines claimed in some bursts? (The 
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evidence for these lines is not very strong, cf. Messina & Share 1992.) (9) Why have 
no galaxies been found in the vicinities of bright 7-ray bursts with well-determined 
positions (Schaefer 1990)? 

The aim of this paper is to describe a qualitative scenario for the production of 
7-ray bursts in mergers of NS-NS and BH-NS binaries, and to demonstrate that there 
are plausible solutions to most of these objections. 

2. MERGER SCENARIO 

The progenitors of close NS-NS and BH-NS binaries must be massive X-ray binaries 
consisting of an O or B main sequence star and a neutron star or a black hole (e.g. Vela 
X-l, Cyg X-l, cf. Trimble 1991 for a review). When the main sequence star evolves, 
the binary very likely undergoes a common envelope phase, after which one has a tight 
binary consisting of the helium core of the OB star and its compact companion. Cyg 
X-3 (Pb=4.79 hr) appears to be an excellent example of this stage of evolution (van 
Kerkwijk et al. 1992). 

We are interested in Cyg X-3 like binaries with separations a ranging from ~ few 
xlO 10 cm, the radius of the helium core (assuming that it behaves like a helium main 
sequence star), to ~ few xlO 11 cm, the limit beyond which gravitational radiation 
losses in the double degenerate binary phase are too slow to cause a merger within 
the age of the universe. To make a NS-NS binary, the helium core needs to have a 
mass Mue between ~ 2.5M© (in order to have a supernova explosion) and 4.2M 
(in order to leave behind a bound NS-NS binary, assuming a symmetric supernova 
explosion and a neutron star mass of 1.4M Q ). After the explosion, we are left behind 
with a NS-NS binary with an orbital eccentricity e = (Mne/2.8) — 0.5, periastron 
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separation a p equal to the pre-explosion separation a Q , and a recoil velocity iWoii = 
180(M H e-2.8)(M H e + 1.4)- 1/2 (ap/10 11 cm)- 1 / 2 km s _1 . Similar estimates for a BH-NS 
binary give somewhat lower orbital eccentricities and space velocities. 

For typical numbers, the recoil velocity of the binary is large enough for the system 
to escape from a small galaxy (though probably not from an L* galaxy). There is 
growing evidence recently that the blue extragalactic light is dominated by faint dwarf 
galaxies between 22 m — 24 m (Cowie 1991, and references therein). These galaxies are 
at modest redshifts (z ~ 0.2 — 0.4) and are very numerous (1 and 6 per square arcmin 
at 22 m and 24 m ). If supernovae are strongly correlated with the blue light, then the 
majority of compact binaries are probably born in low-mass galaxies, escape from their 
hosts, and move distances of up to 1 Mpc in a few xlO 9 yr before they merge. The 
lack of obvious host galaxies associated with 7-ray bursts is then not surprising. The 
burst positions are likely to be offset from that of their parent galaxies by up to an 
arc minute, and at this separation the association will be confused by several other 
faint galaxies in the field. 

The typical duration of a 7-ray burst is ~ 10 s, and the time delay between the 
precursor and main pulse is also ~ 10 s (Murakami et al. 1991). For a NS-NS binary 
in a near-circular orbit, the time to merge is given by t me rge = 38 (a/10 7 5 cm) 4 s = 
160 (Pb/0-ls) 8//3 s. Although the efficiency of tidal interactions and spin-up of the two 
stars as they spiral towards each other is uncertain (Haensel, Paczynski & Amster- 
damski 1991, Kochanek 1992), it is likely that there will be sufficient heating produced 
by tidal dissipation to lead to a super-Eddington luminosity. This may produce some 
of the precursors. The luminosity may also drive a wind from the surfaces of the two 
stars, creating slow-moving baryonic material which may block 7-ray emission from 
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the compact binary. On the other hand, a later-produced relativistically expanding 
burst may be able to shock on the slow wind to produce non-thermal radiation. Yet 
another possibility is that the precursors may in fact be due to the early phases of the 
relativistic burst material plowing through the pre-burst wind. 

Once the two stars merge, there will be ~ few xlO 53 erg of energy divided more- 
or-less equally into three forms: (1) thermal energy, (2) ordered rotational energy, and 
(3) "non-uniform" kinetic energy (either chaotic motions or differential rotation). The 
third form of energy may be quite significant, and may be located quite far from the 
center, if there is strong splashback during the merger as may happen if the rotations 
of the two stars are not synchronized with the orbit. Each of the three energy forms 
will be released through a different channel. 

The thermal energy is radiated mostly as neutrinos and antineutrinos. Goodman 
et al. (1987) showed that in Type II supernovae the finite cross-section associated with 
the reaction, u+u — > e + +e~, results in ~ 10 51 erg being converted into electromagnetic 
energy. Eichler et al. (1989) proposed that this mechanism may produce 7-ray bursts 
in binary mergers. The timescale of the process is the neutrino cooling time, ~ few 
seconds. 

Although modest variability may be expected in neutrino emission because of the 
chaotic (possibly convective) flows in the post-merger fluid, the complicated and highly 
variable time structure observed in some bursts may be too extreme to be explained by 
the neutrino process alone. We suggest that the neutrino mechanism produces the sub- 
class of bursts with relatively smooth profiles and time durations ~ few seconds (e.g. 
GB830801b, Fig. 2 in Higdon & Lingenfelter 1990), and possibly also some precursors. 

The rapid rotation of the post-merger object will almost certainly cause the fluid 
to be dynamically unstable, leading to rapid loss of energy and angular momentum 
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on a timescale of a few msec. This instability dies down when the ratio of kinetic to 
gravitational energy falls below ~ 0.27. At this stage the rotating system will collapse 
directly to a black hole if the total mass is greater than the maximum mass of a neutron 
star (note that some stiff equations of state permit neutron stars more massive than 
2.8©), and if the angular momentum J is smaller than GM 2 /c. If Jc/GM 2 > 1, 
collapse is temporarily prevented, but the system will continue to exhibit a secular 
instability through which it will lose angular momentum on a timescale of a few seconds 
(Friedman 1983). 

Most of the energy released during the first dynamical stage will be in the form of 
gravitational radiation, which could be of importance to future gravitational radiation 
detectors, but we do not expect any direct electromagnetic signal. During the second 
secular phase there may be coherent emission of gravitational radiation that could 
be detected by LIGO, but other processes like neutrino cooling and build up of the 
magnetic field will also take place. 

A fraction of the original mass will not participate directly in the central insta- 
bility and collapse, but will form a surrounding disk/torus which may itself undergo 
dynamical instabilities (e.g. Narayan & Goodman 1989). We consider it likely that 
after this phase a significant amount of residual kinetic energy, < 10 52 erg, will survive 
in a relatively long-lived differentially-rotating disk, and we believe that this energy 
can be tapped almost entirely through electromagnetic processes. 

In the merger of a BH-NS binary, the neutron star will be tidally disrupted, with 
some of the matter being swallowed by the black hole, and the remainder forming a 
torus, generating heat and neutrinos in the process. Subsequent evolution may proceed 
in a similar fashion in both the BH-NS and NS-NS scenarios, and it is not clear at 
this time which of the two is more efficient as a 7-ray burster. 
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The mechanism we propose to produce 7-ray bursts similar to GB820331 and 
GB841215 (Fig. 2 in Higdon & Lingfelter 1990), which have complex time profiles, is 
based on magnetic fields. Immediately after the merger the field strength is probably 
~ 10 12 G, but instabilities such as the Balbus & Hawley (1991) mode will cause the 
field to build up quickly to something like an equipartition value ~ 10 16 — 10 17 G. 
Even if there is no instability, the shearing action of the differential rotation will lead 
to field build-up on a slower timescale ~ 10 4 rotations. (Note that, while electric 
fields are limited to a maximum strength of ~ 10 13 V cm -1 by the pair creation 
threshold, magnetic fields can build up to much higher strengths.) Once the field 
achieves equipartition, it will tend to exhibit the Parker (1966) instability in which 
the field will float up and break out of the disk (relativistically in our case) on a 
dynamical time, leaving the matter behind. Particularly near the top of the magnetic 
loop, the baryonic contamination may be quite low, making the conditions favorable 
for a gamma-ray burst. The burst itself will probably be created as a series of explosive 
reconnection events in the rising magnetic field, just as in solar flares. For the assumed 
field strength, the temperature of the flare will be ~ few xlO 11 K, the photons will 
have characteristic energies ~ 10 MeV, and the total energy in a single burst will be 
~ 10 50 erg for linear dimensions ~ 1 — 10 km. The resulting optically thick fireball 
will expand relativistically as described by Paczyhski (1986), Goodman (1986), Shemi 
& Piran (1990), and Paczyhski (1990). 

A very important feature of the model is that the rate of energy release is not 
limited by the Eddington luminosity, but by the efficiency with which the Parker 
instability can separate magnetic fields from matter. The energy dissipation occurs 
above the disk, and, because of the relativistic outward motion of the flare, very little 
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energy returns to the disk. If this aspect of the scenario works efficiently (which needs 
to be demonstrated), then we have a viable solution to the problem of how to separate 
the photons from the baryons (how to "separate light from the darkness"). 

The flare activity will make the burst profile quite complex, with many sub- 
bursts. Also, the total duration over which the successive sub-bursts occur will be 
quite variable from one object to another, depending on the details of the post-merger 
fluid configuration. These features eliminate some of the strongest objections to the 
merger scenario. The important point is that the model has two distinctly different 
timescales: firstly, a dynamical time, ~ 10~ 3 s, which may be related to the rise 
times of individual flares (because the Parker instability, once it gets going, operates 
essentially on a dynamical time), and, secondly, an accretion (or magnetic "viscosity") 
time, ~ 1 — 100 s, related to the duration of the whole 7-ray burst. We cannot calculate 
the second timescale from fundamental theory, but we note that in other objects with 
disks, such as cataclysmic variables and X-ray binaries, the accretion time is typically 
many orders of magnitude longer than the dynamical time. 

A generic problem with any cosmological model is that the huge initial optical 
depth leads to themalization of the fireball (Paczynski 1990) and a black body spec- 
trum. This is in clear conflict with the observed spectra of 7-ray bursts, which have a 
significant excess of X-rays and very hard 7-rays compared to any Planck spectrum. 
There are at least two possible solutions to this problem, both related to the flare- 
like energy release and rapid variability. First, while the instanteneous spectrum of a 
local flare may be Planckian, the superposition of many Planck curves with a wide 
range of temperatures and intensities may appear as a very broad broken power law. 
Second, it is likely that the flares will eject some mass due to imperfect separation of 
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energy and matter. These relativistic ejecta may collide with each other, as well as 
with the pre-merger wind, at a fairly large radius and at low optical depths, giving 
a non-Planckian spectrum through various non-thermal processes. For instance, given 
the strong magnetic fields, synchrotron processes might naturally produce the observed 
power-law spectrum if enough relativistic particles are present. 

An important point to note is that the observed photon energy /iz/ Q b s * s related to 
the emitted energy in the frame of the radiating fireball hv era by the large relativistic 
7 ~ 10 2 — 10 3 of the expansion (Paczynski 1986, Goodman 1986). Therefore, there is 
no difficulty in producing very hard spectra. In particular, note that a cut-off at 511 
keV in the emitter frame translates to a cut-off at 7 x 511 keV for the observer, which 
can correspond to several hundred MeV. 

The fact that the energy release in the model is predominantly through a magnetic 
process may solve another potential problem, the cyclotron absorption lines that have 
been claimed in some bursts. In very rough terms, an initial field of 10 16 G at the 
point of initiation of the burst will be diluted to a field strength ~ 10 10 G in the frame 
of the emitting fluid if the expansion factor is ~ 10 3 . However, because of relativistic 
beaming, the observed cyclotron line will appear as if it is emitted in a field ~ 7 x 10 10 
G ~ 10 13 G, which is roughly consistent with the claimed field strengths (Murakami et 
al. 1988, Fenimore et al. 1988). The narrowness of the observed lines requires a stable 
value of the Lorentz 7 for the duration of the line. This will be difficult to maintain 
for any length of time in our highly erratic burster, in qualitative accordance with the 
"now you see them, now you don't" nature of the reported lines. 

Another interesting possibility should be finally mentioned. As the disk cools 
through neutrino emission, it might at some point make a transition to a supercon- 
ducting superfluid state. Plasma processes in a degenerate superconducting superfluid 

10 



have been hardly studied at all. It is conceivable that the separation of the field and 
the baryons may be particularly efficient in a magnetic superfluid. This problem merits 
further investigation. 

3. CONCLUSIONS 

We propose a scenario for gamma-ray bursts involving merging NS-NS and BH-NS 
binaries at cosmological distances. In our picture, complex main bursts with highly 
variable time profiles are produced through many successive magnetic flares and re- 
connection events tapping the free energy in differential rotation in a post-merger disk. 
Less variable bursts may be produced by neutrinos converting a fraction of their en- 
ergy into electromagnetic radiation. The precursors of the bursts may be produced by 
the neutrino mechanism, or by the initial phases of the post-merger ultra-relativistic 
burst ejecta plowing through the slow wind created in the pre-merger tidal heating 
phase, or even possibly by the pre-merger heating itself. This scenario is capable of 
explaining the qualitative features of most of the observations, and circumvents many 
of the arguments against the cosmological merger model of gamma-ray bursts. 

The model employs two kinds of objects, NS-NS and BH-NS binaries, and two 
distinct mechanisms, magnetic flares and neutrino interactions. Different combinations 
of these may have distinct observational signatures which could be potentially identified 
in the large database of burst profiles being acquired by BATSE. 

Apart from 7-rays, it is likely that the bursts in our scenario will also eject ~ 
10 51 erg of energy in ultra-relativistic ejecta, i.e. cosmic rays, with little or no non- 
relativistic matter. These ejecta should create a supernova-like remnant far from the 
parent galaxy, where the density of diffuse matter may be very low. 
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The scenario makes a definite prediction. If and when LIGO is commissioned (Vogt 
1992), strong 7-ray bursts should be accompanied by gravity wave detections (though 
the reverse need not necessarily be true, particularly if the 7-rays are beamed). Since 
the strength of the gravitational radiation signal from a pre-merger binary can be 
calculated accurately (Lincoln & Will 1990) except for the unknown inclination, LIGO 
should provide good distance estimates to individual bursts and also pinpoint the exact 
time of the merger. This information will permit more detailed interpretation of the 
7-ray burst observations. 

The scenario we have presented in this paper is at this point no more than a 
qualitative sketch of some likely possibilities. The purpose here is merely to show that 
none of the many arguments against the merger model are necessarily fatal. More 
detailed analyses of the individual ideas are needed now. 

One other point worth emphasizing is that the scenario is in a real sense the 
most conservative one possible at the present time. Based purely on the observations, 
a cosmological location for 7-ray bursts is the most reasonable hypothesis; all other 
proposals need coincidences at some level in the properties of the associated source 
population. Within the cosmological framework, our scenario employs the most obvious 
source population that we can think of, double compact binaries. This population is 
definitely known to exist, we know its members will merge, and we can be certain that 
huge quantities of energy will be released. However, even if this particular model turns 
out to be incorrect, the cosmological hypothesis itself will survive, given the observed 
distribution of burst positions and intensities. 

We thank Roger Blandford and John Friedman for discussions and comments. 
This work was supported by NASA grants NAG 5-1901 and NAG 5-1904. 
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